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ABSTRACT 

This  report  presents  an  analytical  method  for  computing 
the  forces  and  moments  on  a  rigid  foil  of  finite  span  undergoing 
large  amplitude  pitching,  heaving,  and  surging  oscillations.  The 
influence  of  foil  stall  and  the  induced  slipstream  are  included. 
Experimental  data  obtained  on  a  large  amplitude  oscillating  foil 
propulsor  are  presented  and  compared  with  the  theory.  Perfor¬ 
mance  predictions  for  an  oscillating  foil  propulsor  for  use  on  a 
small,  15  knot,  shallow-draft  boat  are  also  presented. 

INTRODUCTION 

Propulsion  of  a  vehicle  in  weedy  shallow  waters  by  means  of 
a  rotating  propeller,  even  when  surrounded  by  a  nozzle  or  enclosed 
within  a  casing,  as  in  a  pump-jet  system,  is  very  difficult  as 
weeds  tend  to  block  the  inlet  to  the  system  and  wind  around  the 
rotating  shaft  of  the  impellers.  A  possible  scheme  for  alleviating 
the  problem  is  to  achieve  propulsion  with  an  oscillating  foil. 
Entanglement  will  be  avoided  since  no  continuously  rotating  shaft 
is  exposed  to  the  weeds. 

In  order  to  achieve  practical  levels  of  thrust  such  an  oscil¬ 
lating  foil  must  undergo  large  amplitude  oscillations  at  rela¬ 
tively  high  frequency.  The  magnitude  of  these  oscillations  is 
such  that  the  classical  small  amplitude  theories  cannot  provide 
adequate  performance  predictions  for  engineering  purposes. 
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This  study  presents  an  analytical  method  for  computing  the 
forces  and  moments  on  a  rigid  foil  of  finite  span  undergoing  large 
amplitude  pitching,  heaving,  and  surging  oscillations.  The  method 
is  capable  of  dealing  with  cambered  foils  and  with  motions  where 
the  mean  angle  of  attack  is  not  zero.  The  influence  of  foil  stall 
and  the  induced  slipstream  has  also  been  included. 

An  experimental  program  was  undertaken  to  determine  the  va¬ 
lidity  of  the  analysis.  It  was  found  that  good  agreement  was  ob¬ 
tained  under  most  of  the  test  conditions  except  for  the  case  of 
zero  forward  speed.  Under  this  condition  the  predicted  thrust 
was  less  than  that  measured. 

The  analytical  procedure  was  used  for  the  preliminary  design 
of  an  oscillating  foil  propulsor  suitable  for  use  on  a  15  knot, 
shallow-draft  boat  of  2000  pound  payload,  such  as  the  "ski-barge." 
The  results  indicate  that  an  oscillating  foil  propulsor  can  provid 
efficient  shallow  water  propulsion  with  a  high  degree  of  maneuvera 
bility.  The  ultimate  practicality  of  the  system  will  depend  more 
on  the  solution  of  the  mechanical  drive  problems  than  with  the 
hydrodynamic  performance. 

This  study  of  oscillating  foil  propulsors  was  carried  out  for 
the  U.  S.  Army  Engineering  Research  and  Development  Laboratories, 
Ft.  Belvolr,  Virginia,  under  Contract  Number  DA-44-009-AMC-1759(T ) 
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THEORY 

A  method  for  estimating  the  forces  and  moments  on  a  foil 
performing  large  amplitude  oscillations  is  given  in  the  following 
section.  For  convenience  the  equations  have  been  summarized  in 
Appendix  I.  The  equation  numbers  appearing  here  are  thus  those 
of  Appendix  I. 

The  problem  considered  is  that  of  a  rigid,  finite  span  foil 
supported  in  a  uniform  stream  and  undergoing  regular  pitching, 
heaving,  and  surging  oscillations.  It  is  not  required  that  the 
foil  be  symmetric  or  that  the  mean  angle  of  attack  be  zero.  The 
phase  relation  between  the  pitching,  heaving,  and  surging  motions 
is  arbitrary,  but  they  must  all  be  of  the  same  frequency.  It  is, 
however,  necessary  to  know  the  static  performance  of  the  foil. 

The  motions  and  forces  are  defined  at  an  arbitrary  point  on  the 
chord  (or  its  extension)  of  the  mean  aerodynamic  foil  section. 

This  results  in  no  loss  of  generality  because  the  inclusion  of 
surging  oscillations  makes  it  possible,  through  a  transfer  of  co¬ 
ordinates,  to  compute  the  forces  and  moments  on  a  foil  whose  pitch 
axis  lies  above  or  below  the  foil  .chord. 

Method  of  Analysis 

The  hydrodynamic  analysis  is  based  on  the  classical  small 
amplitude  theory  of  Theodorsen  (Reference  6)  and  Garrick  (Refer¬ 
ence  3)  for  a  foil  oscillating  with  sinusoidal  motions  in  a  free 
stream.  This  theory  was  first  divided  into  terms  which  could  be 
associated  with  "added  mass"  and  hydrodynamic  circulation.  The 
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important  coefficients  were  then  Identified  so  that  their  theo¬ 
retical  linear  values  could  be  replaced  with  more  exact  values. 

It  was  found  that  the  forces  and  moments  due  to  hydrodynamic 
circulation  at  any  instant  of  time  could  be  considered  to  be  equal 
to  the  steady-state  values  that  would  occur  at  an  equivalent  angle 
of  attack  ag  and  velocity  V,  plus  an  additional  force  and  moment 
which  result  from  the  angular  rotation  of  the  foil.  This  is  for¬ 
tunate  since  it  makes  it  possible  to  utilize  the  wealth  of  existing 
three-dimensional  foil  data  and  theory.  Finally,  the  linearized 
values  of  the  characteristic  angles,  velocities,  and  acce] era tions 
were  replaced  by  their  actual  values.  The  inertial  forces  arising 
from  the  foil  mass  are,  of  course,  obtained  from  rigid  body  me¬ 
chanics  . 

The  resulting  equations  are  summarized  in  Appendix  I  while 
the  derivation  of  the  equations  for  the  hydrodynamic  force  is 
given  in  Appendix  II. 

The  equations  reduce  to  those  of  the  classical  linear  theory 
for  all  frequencies  when  the  amplitude  of  the  motions  is  small  and 
there  is  no  surge.  They  also  reduce  to  the  quasi-steady-state  so¬ 
lution  for  all  amplitudes  when  the  frequency  is  low.  While  surg¬ 
ing  motions  are  not  included  in  the  original  theory,  they  appear 
here  simply  as  a  result  of  permitting  large  amplitude  motions. 

The  important  assumptions  underlying  the  analysis  are  that 
the  vorticity  is  distributed  sinusoidally  in  the  wake  and  that  the 
wake  lies  in  an  infinitely  thin  plane  containing  the  foil  chord. 
Obviously,  the  wake  will  not  lie  in  a  thin  plane  but  will  be  dis¬ 
tributed  vertically  through  a  height  equal  to  the  vertical  foil 
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displacement.  Further,  the  nonlinear! ties  introduced  by  permitting 
large  amplitude  motions  will  result  in  a  wake  that  does  not  have 
a  sinusoidal  vortex  distribution  even  when  the  foil  motions  are 
sinusoidal.  However,  it  is  felt  that  this  will  not  Introduce  ap¬ 
preciable  errors  as  long  as  either  the  wavelength  of  the  vortex 
wake  is  large  compared  to  the  foil  chord  (low  frequency)  or  the 
amplitude  of  oscillations  is  small  compared  to  the  foil  chord.  By 
the  same  token,  little  error  should  be  introduced  if  the  motions 
are  not  exactly  sinusoidal.  Errors  will,  of  course,  be  introduced 
for  foil  motions  which  are  of  both  large  amplitude  and  high  fre¬ 
quency.  The  experiments,  discussed  later  in  this  report,  indi¬ 
cate  that  under  these  conditions  the  theory  tends  to  underesti¬ 
mate  the  forces  on  the  foil. 

A  further,  arbitrary,  assumption  concerns  the  unsteady  stall¬ 
ing  characteristics  of  the  foil.  It  has  been  assumed  that  the 
foil  will  stall  at  a  lift  coefficient  yielding  the  same  strength 
leading  edge  vorticity  that  would  occur  at  the  steady-state  stall 
lift  coefficient.  This  is  equivalent  to  assuming  that  leading 
edge  stall  occurs.  This,  however,  is  not  an  inherent  limitation 
and  can  be  revised  when  more  experimental  data  is  available.  The 
nature  of  the  stall  will  obviously  depend  on  the  foil  shape,  sur¬ 
face  condition,  and  instantaneous  Reynolds  number.  It  would  in¬ 
deed  be  fortunate  if  such  a  simple  relation  existed  between  the 
steady-state  and  unsteady  stall  lift  coefficients.  The  method 
of  calculating  the  unsteady  stall  lift  coefficient  will  be  dis¬ 
cussed  later. 
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Although  the  present  analysis  considers  only  a  rigid  foil, 
the  method  can  easily  be  extended  to  include  a  foil  with  a  flap 
since  this  problem  has  already  been  investigated  by  Garrick 
(Reference  3 ) • 

It  must  be  made  clear  that  the  present  analysis  is  not  a 
true  nonlinear  theory.  It  is,  rather,  a  combination  of  the  non- 
linearities  obtained  from  a  quasi-steady  state  analysis  with  the 
results  of  a  small  amplitude  unsteady  analysis.  As  such,  it  is 
subject  to  certain  limitations  imposed  by  the  assumptions  involved 

Forces  and  Moments  on  an  Oscillating  Foil 

Consider  a  foil  with  a  mean  aerodynamic  chord,  c,  supported 
in  a  stream  of  velocity  Uq.  The  foil  is  performing  simple  har¬ 
monic  pitching,  heaving,  and  surging  oscillations  at  a  rate  of  o> 
radians  per  second  with  a  mean  angle  of  attack  am>  These  motions 
are  defined  at  a  point  0  located  a  distance  d  ahead  of  the  quarter 
chord  of  the  mean  aerodynamic  chord.  Figure  1  illustrates  the 
geometry  involved.  We  wish  to  determine  the  lift,  drag,  and  mo¬ 
ment  at  0  as  a  function  of  time.  At  any  instant  of  time  t,  it 
is  convenient  to  separate  these  forces  into  contributions  from 
hydrodynamic  circulation,  added  mass,  and  real  foil  mass.  Thus, 
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where  the  subscripts  "c" ,  "am",  and  "m"  refer  to  the  contribu¬ 
tions  from  circulations,  added  mass,  and  foil  mass  respectively. 
The  first  two  contributions  are  obtained  from  the  hydrodynamic 
analysis  and  will  be  dealt  with  first.  The  last  contribution 
is  just  the  inertial  forces  arising  from  the  foil  mass  and  will 
be  dealt  with  later. 

The  oscillating  lift  on  the  foil  results  in  the  shedding 
of  an  oscillating  vortex  wake.  This  wake  induces  a  flow  at  the 
foil  which  is  proportional  to  the  foil  lift  and  a  function  of  the 
number  of  chord  lengths  the  foil  has  traveled  through  the  fluid 
in  one  cycle.  This  distance  can  be  expressed  as  ir/ k  where  k  is 
the  reduced  frequency,  and  is  defined  as 


k 


O)  c/2 
U 


[8] 


Here  U  is  the  total  streamwise  velocity  at  the  foil  and  includes 
the  induced  velocity  in  the  streamwise  direction  u.  Thus, 


U  =  U  +  u 
o 

Since  it  is  not  possible  to  compute  the  slipstream  velocity  at 
the  foil  until  the  foil  loading  is  known,  the  total  velocity  U 
is  assumed  known.  The  mean  induced  velocity  u  is  then  computed 
and  the  free  stream  velocity  is  taken  to  be 

UQ  =  U  -  u  [29] 

The  method  of  computing  u  will  be  discussed  later. 


f 
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The  induced  flow  is  characterized  by  the  complex  function 


C(k) 


C(k)  =  P(k)  +  i  G (k ) 


[7] 


Values  of  P  and  G  were  computed  by  Theodorsen  for  infinite  span 
and  have  since  been  computed  by  others  for  finite  span  foils  and 
foils  under  a  free  surface  (References  1,  2,  and  4).  Values  of 
F  and  G  are  plotted  in  Figure  2A  for  aspect  ratios  from  0  to  10. 
(These  plotted  values  have  been  normalized  by  dividing  by  their 
steady  state  values  so  they  must  be  used  in  conjunction  with  the 
corresponding  three-dimensional,  steady  state  foil  performance). 

Because  the  flow  induced  by  the  wake  is  cyclic  there  is  at 
any  Instant,  both  an  induced  velocity  and  an  induced  acceleration 
relative  to  the  foil.  In  general,  the  induced  velocity  will 
modify  the  instantaneous  angle  of  attack  while  the  induced  ac¬ 
celeration  will  modify  the  inertial  forces  from  "added  mass". 
Since  the  strength  of  the  wake  vorticity  is  proportional  to  the 
strength  of  the  oscillating  circulation  on  the  foil,  there  is  a 
coupling  between  the  inertial  and  circulating  forces.  It  is 
mathematically  equivalent,  however,  to  treat  the  influence  of  the 
wake  as  modifying  the  angle  of  attack  and  causing  a  time  lag  in 
the  circulatory  forces.  This  is  a  useful  interpretation  since 
it  uncouples  the  inertial  forces  from  the  circulatory  forces 
and  permits  the  use  of  an  "equivalent  steady  state  angle  of 
attack"  for  computing  the  circulatory  forces.  These  notions  are 
important  in  the  present  analysis  because  they  permit  the  identi¬ 
fication  and  separation  of  many  important  coefficients  so  that 


I 


'  iJ 
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their  linearized  values  can  be  replaced  by  nonlinear  or  experi¬ 
mentally  obtained  values.  The  uncoupling  of  the  linearized  equa¬ 
tions  and  the  identification  of  the  various  coefficients  is 
presented  in  Appendix  II. 

The  unsteady  forces  and  moments  due  to  hydrodynamic  circula¬ 
tion  at  any  time  t  are  equal  to  the  steady-state  values  that  would 
occur  at  an  equivalent  angle  of  attack  and  velocity  V,  plus  an 
additional  force  and  moment  which  result  from  the  angular  rotation 
of  the  foil.  V  is  taken  as  the  instantaneous  velocity  at  the  mid¬ 
chord  and  is  given  by  Equation  [6]  in  Appendix  I.  ag  is  composed 
of  the  mean  angle  of  attack  and  an  equivalent  unsteady  angle 

of  attack  a  1 .  a  '  is  numerically  equal  to  the  oscillatory  angle 
0  0 

of  attack  at  the  3/^  chord  evaluated  at  time  (t  -  At)  and  multi¬ 
plied  by  the  factor  "V  F2  +  G2 .  Thus, 

a  -  a  +  a  1 
e  m  e 

where 

a  '  =  a  /,  ~\/ F2  +  G2 
e  3  c/4  V 


[10] 


[11] 


a3c/4  lns't:an^aneous  oscillatory  angle  of  attack  at  the 

3/4  chord  and  is  given  by  Equation  [12]  in  Appendix  I.  The  bar 
is  used  to  signify  that  this  quantity  is  evaluated  at  time  (t-At) 
where 


At  = 


arctan  (-G/F) 

CD 


[14] 
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The  effect  of  the  vortex  wake  on  the  circulatory  forces 
associated  with  angle  of  attack  can  therefore  be  interpreted  as 
multiplying  the  oscillatory  portion  of  the  angle  of  attack  by 
V P2  +  G2  and  introducing  a  phase  lag  of  oAt  with  respect  to  the 
angle  of  attack  at  the  3/^  chord.  Thus  the  oscillatory  portion 
of  the  vortex  wake  introduces  an  induced  angle  of  attack  of 


a3/4  c  [1-0(k>]  "  a3/it  o  '  ae ' 

Values  of  V  F2  +  G2  and  o>At  are  given  in  Figure  2B  for  aspect 
ratios  from  0  to  10.  The  average  angle  of  attack  is  unaffected. 
The  equivalent  steady-state  forces  and  moments  at  the  1/4  chord 
can  now  be  found  from  the  following  relations 


\ 

Li  =  CL  ipV2  S 


Di 


cD  ipV2  S 
CM  |pV2Sc 

Min 

h  ° 


I 


[3] 


These  are  to  be  evaluated  from  the  static  foil  performance  at 
a  =  ag.  Li  and  Di  are  found  to  act  normal  and  parallel  to  the 
instantaneous  stream  direction  at  the  quarter  chord  just  as  in 
the  case  of  steady  flow.  If  the  foil  has  angular  rotation  the 
angle  of  attack  will  vary  along  the  chord  by  an  amount  y  per 
quarter  chord 
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7 


a  r  c  ta  n 


c  a 
4  V 


[16] 


The  angle  of  attack  at  the  1/4  chord  Aa  can  then  be  expressed  as 


Aa  =  a  -  2y  [ 15 ] 

e 

where  2y  Is,  of  course,  just  the  difference  in  angle  of  attack 
between  the  3/4  chord  and  the  1/4  chord. 

The  presence  of  foil  rotation  also  gives  rise  to  an  addi¬ 
tional  force  Si  and  moment  M  .  Si  is  a  thrust  force  attributed 

7 

to  "leading  edge  suction"  caused  by  additional  vorticity  at  the 
leading  edge.  This  vorticity  is  due  to  an  increase  in  velocity 
around  the  leading  edge  which  results  from  the  angular  rotation. 
Si  acts  in  the  plane  of  the  chord  and  can  be  expressed  as 


Si  =  2iry  ipV^S  tan  y  [4] 

M  is  a  damping  moment  in  that  it  always  opposes  the  angular  ro¬ 
tation.  It  can  be  expressed  as 

My  =  -  |  7  ipV2Sc  [51 

The  terms  Si  and  M  are  of  special  interest  since  they  do 

7 

not  appear  when  there  is  no  foil  rotation.  These  terms  are  un¬ 
affected  by  the  vortex  wake  and  are  thus  independent  of  the  foil 
lift,  the  past  history  of  the  motion,  and  the  aspect  ratio.  An 
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energy  balance  reveals  that 
foil  angular  rotation  (M  &) 
by  the  thrust  force  (SiV). 


the  work  required  to  maintain  the 
is  totally  recovered  in  the  work  done 
That  is 


M  a  =  SiV  =  -  y  4pV2Scd 
Y  2 

This  is  an  important  result  for  oscillating  foil  propulsion  for 
two  reasons.  First  of  all,  the  propulsive  efficiency  from  this 
source  at  any  instant  is  equal  to  U^/Q  cos  a.  Since  a  reaches 
its  peak  value  near  the  end  of  the  stroke  when  a  is  small  (or 
zero)  the  propulsive  efficiency  from  this  source  will  be  very 
high.  Secondly,  the  propulsive  force  due  to  Li  goes  to  zero  at 
the  end  of  each  stroke  but  again  this  is  when  Si  reaches  its  maxi¬ 
mum.  Thus  the  presence  of  thrust  from  Si  will  tend  to  smooth  out 
the  thrust  variations  and  increase  the  overall  thrust  coefficient. 

Of  course,  if  y  becomes  too  large  the  flow  will  separate 
from  the  leading  edge  and  the  foil  will  stall.  In  general  the 
nature  of  the  stall  will  depend  on  the  foil  shape,  surface  con¬ 
dition,  instantaneous  Reynolds  number,  and  frequency  and  magnitude 
of  tiie  oscillation.  However,  if  it  is  assumed  that  this  leading 
edge  stall  will  occur  when  the  vorticity  at  the  leading  edge 
reaches  the  same  value  as  when  steady-state  stall  occurs,  then 
the  steady  state  stall  lift  coefficient  can  be  used  to  estimate 
when  stall  will  occur  on  the  oscillating  foil.  If  we  define  a 
leading  edge  suction  coefficient  Cg  as 


^pV2S  tan  y 


[4] 
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then  the  strength  of  the  instantaneous  leading  edge  vorticity 
will  be  the  same  as  for  a  steady  foil  when 


C  ) 
s  os 


cilia  ting 


steady 


If  a  maximum  and  minimum  lift  coefficient  can  be  determined  from 
the  steady  foil  performance  then  the  following  criteria  can  be 
applied  to  the  unsteady  foil  to  estimate  when  stall  will  occur. 


If  (C  -  C  )  <  0 

ij  s  mm 


or  (C  -  C  )  >  C  the  foil  will  stall, 

s  b  max 


Of  course,  if  a  and  a  are  both  of  the  same  sign  than  C 

e  s 

will  tend  to  offset  CT  and  permit  values  of  C  beyond  those  where 

.L  Li 

stall  would  occur  under  steady  conditions.  However,  trailing 

edge  stall  may  eventually  put  an  upper  limit  on  C  .  The  magni- 

L 

tude  of  this  upper  limit  cannot  be  determined  from  steady  tests 

since  it  involves  the  severity  of  the  adverse  pressure  gradient 

at  the  trailing  edge.  It  can  probably  only  be  determined  from 

tests  on  an  oscillating  foil.  It  seems  likely,  however,  that 

instantaneous  maximum  values  of  C  could  be  as  large  as  twice  the 

L 

maximum  steady  value  and  this  has  been  assumed  in  the  absence  of 
experimental  data.  Combining  this  criteria  with  the  previous 
criteria  yields  the  conditions 


If  t(C 


or 


(0/2))  <  C  . 

^  Jj  mm 


or  [  (C 


C  )  or  (CT/S)  ]  >  C 


then. 


max 


I 
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the  foil  will  stall  and  the  calculated  values  of  CT  and  C  are 
invalid.  If  this  occurs  then  some  estimate  must  be  made  of  the 
stalled  values  of  Li,  Do.,  Mi,  Si  and  M  .  The  values  assumed  for 
these  quantities  after  stall  are  discussed  later  in  this  report 
in  conjunction  with  the  computer  solution  of  the  problem.  Con¬ 
sidering  the  complicated  physics  involved  in  determining  when 
stall  will  occur  in  an  unsteady  flow,  it  would  indeed  be  fortunate 
if  such  simple  relations  existed  between  the  steady-state  and  un¬ 
steady  stall  lift  coefficients.  The  above  criteria  should  be  con¬ 
sidered  only  as  an  approximate  method  of  estimating  stall  and 
should  be  revised  when  data  or  a  more  exact  method  are  available. 
The  data  presented  in  this  report  indicate  that  these  criteria 
are  probably  reasonable.  However,  the  nature  and  duration  of  the 
stall  cannot  be  determined  from  these  data  since  only  mean  values 
were  recorded. 

It  is  also  of  interest  to  note  the  similarity  of  these  forces 
and  moments  with  those  which  occur  on  a  steady  foil  with  circular 
arc  camber  of  camber-chord  ratio  equal  to  y/2.  It  should  first 
be  noted  that  the  change  in  angle  between  the  camber  line  and 
the  free  stream  direction  for  the  cambered  foil  is  also  equal  to 
y  per  quarter  chord  as  can  be  seen  in  Figure  3 •  For  such  a  foil 
the  zero  lift  line  is  parallel  to  the  camber  line  at  the  3/^  chord 
so  that  the  effective  angle  of  attack  is  measured  as  the  angle 
between  the  stream  and  the  3/^  chord  just  as  in  the  present  case. 
In  addition  the  1/4  chord  pitching  moment  for  the  cambered  foil 
is  also  given  by  Equation  [5].  The  force  Si,  however,  does  not 
exist  on  the  steady,  cambered  foil.  The  oscillating  foil  thus 


I 
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behaves  much  like  a  cambered  foil  at  an  angle  of  attack  Aa  with 
sufficient  camber  to  generate  the  same  lift  as  a  flat  foil  at  an 
angle  of  attack  (i.e.  the  effective  angle  of  attack  due  to  foil 
rotation  is  2y) . 

The  forces  and  moments  from  circulation  can  now  be  transferred 

from  the  quarter  chord  to  the  reference  point  0  and  resolved  so 

that  the  lift  L  is  normal  to  the  free  stream,  and  drag  D  is  paral- 

o  c 

lei  to  the  free  stream.  This  transformation  is  accomplished  by 
the  following  equations 


Lc  =  Li  cos  (a-Aa)  -  Di  sin  (a-Aa)  4-  Si  sin  (a) 

D  =  Di  cos  (a-Aa)  +  Li  sin  (a-Aa)  -  Si  cos  (a)  >  [2] 

c 

=  Mi  -  d  [ Li  cos  (Aa)  +  Di  sin  (Aa)]  +  M 

Forces  and  moments  associated  with  added  mass  are  unaffected 
by  the  wake  and  are  simply  proportional  to  the  product  of  the  foil 
acceleration  relative  to  the  fluid  and  the  added  mass  acting  in 
the  direction  of  the  acceleration.  These  forces  and  moments  are 
given  by  the  following  relations 


Ijp  =  atn.w. 


D2  =  am  » ( t/c  )  w 


[19] 


M2  =  -am'K(c/4)  a 
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La  and  D2  are  found  to  act  through  the  midchord  normal  and  paral¬ 
lel  to  the  chord  line  respectively.  M2  is  the  midchord  pitching 
moment.  The  term,  "am",  is  the  added  mass  of  a  foil  heaving  nor¬ 
mal  to  the  plane  of  its  chord.  Reference  5  gives  an  approximate 
formula  for  determining  "am"  for  finite  aspect  ratio  wings  of 

am  =  — == .  p7r(c/2)2b  [20] 

Vl+A2 

The  effective  added  mass  acting  in  the  chordwise  direction  is  as¬ 
sumed  to  be  (t/c)2  times  the  added  mass  acting  normal  to  the  chord 

The  quantities  and  wc  are  the  components,  normal  and  paral 
lei  to  the  chord  respectively,  of  the  rate  of  change  of  velocity 
relative  to  the  midchord.  These  are  given  by  Equation  [21]  in  Ap¬ 
pendix  I.  The  product  [am  K(c/4)2]  in  the  equation  for  M2 
(Equation  [19])  Is  the  moment  of  inertia  of  the  added  mass.  The 
effective  radius  of  gyration  of  the  added  mass  is  therefore 
V"k(c/4).  Typically,  K  has  values  of 


Approximate  values  of  K  can  be  found  for  other  aspect  ratios  with 
the  following  equation 


0. 17A 
~  A+1.43 


+  0.33 


[  22B] 


J 


4 


The  forces  and  moments  from  added  mass  can  be  transferred 
from  the  midchord  to  the  reference  point  0  by  the  following  re¬ 
lations 

L  =  L2  cos  a  -  D2  sin  a 
am 

D  =  D2  cos  a  +  L2  sin  a 
am 

M  =  M2  -  (c/4  +  d)L2 

Inertial  forces  due  to  the  actual  foil  mass  are  found  from 
rigid  body  mechanics  and  are  given  at  the  reference  point  0  by 
the  following  equations 


Lm  ~  _m[h  +  (c/4  +  d-s)(a2sin  a  -  a  cos  a)] 

•• 

Dm  =  m[^  +  (c/4  +  d-s)  (a2cos  a  +  a  sin  a)] 

M  =  -  m  r2  a  -  (c/4  +  d-s)(L  cos  a  +  D  sin  a) 
m  '  /  '  m  m 


[23] 


Here  m  is  the  actual  foil  mass,  s  is  the  location  of  the  center 
of  gravity  ahead  of  the  foil  midchord,  and  r  Is  the  radius  of 
gyration  about  the  center  of  gravity. 

There  is  a  basic  difference  between  the  inertial  forces  as¬ 
sociated  with  added  mass  and  those  associated  with  real  mass.  The 
inertial  forces  from  added  mass  are  proportional  to  the  foil  ac¬ 
celeration  relative  to  the  fluid  while  those  from  real  mass  are 
proportional  to  the  foil  acceleration  relative  to  an  inertial 
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coordinate  system.  As  a  result,  the  added  mass  experiences  no 
centrifugal  force  as  in  the  case  of  the  real  mass.  (The  centri¬ 
fugal  force  appears  as  the  (a)2  terms  in  Equation  [23]). 

The  mean  values  of  the  various  forces  and  moments  are,  of 
course,  obtained  by  integrating  these  forces  and  moments  over  an 
integral  number  of  cycles  and  dividing  by  the  total  time  of  in¬ 
tegration.  Thus 


L  = 


2n7r 


2n7T 


D  = 


2nm 


D  d(oyfc) 


2n7r 


M 


2  mr 


M  d  (cot ) 


o 


where  n  is  the  number  of  cycles. 

In  the  case  of  a  propulsion  system  we  are  interested  in  the 
mean  thrust  T  and  mean  power  p.  Thrust,  of  course  is  just  the 
negative  of  drag  while  the  power  is  made  up  of  the  rate  at  which 
L,  D,  and  M  do  work 


[26] 


P  =  - 


2nn 


2nir 

(Lh  -  D£  +  Ma  )d  (cot ) 


[27] 


o 


HYDRONAUTICS,  Incorporated 


-19- 


Slipstream 

It  Is  to  be  expected,  from  the  basic  principle  of  action 
and  reaction,  that  when  the  foil  has  a  net  positive  thrust  there 
must  be  a  reaction  on  the  fluid  which  pushes  it  backward  and  forms 
a  slipstream.  This  mechanj sm  has  been  described  by  von  Karman  and 
Burgers  on  page  308  of  Reference  7.  When  the  foil  is  producing 
thrust,  the  circulation  is  greater  on  the  downward  stroke  than  on 
the  upward  stroke.  (Positive  circulation  being  that  which  pro¬ 
duces  an  upward  lift).  Hence  when  the  foil  is  at  its  highest 
point,  the  circulation  will  be  increasing,  and  thus  a  vortex  will 
detach  itself  from  the  trailing  edge  with  rotation  opposite  to  the 
direction  of  the  increasing  circulation.  When  the  foil  is  at  its 
lowest  position,  the  circulation  is  decreasing  and  a  vortex  will 
detach  with  the  opposite  rotation.  The  vortex  sheet  thus  shed  has 
the  form  of  a  wave  with  a  wavelength  of  v/k  chords  in  which  the 
vorticity  varies  across  the  wake  from  a  negative  strength  at  the 
top  to  a  positive  strength  at  the  bottom.  This  is  illustrated  in 
Figure  4.  Because  of  the  wavey  nature  of  this  vortex  sheet  an  in¬ 
duced  velocity  in  the  downstream  direction  is  generated.  This  in¬ 
duced  velocity  has  a  form  similar  to  a  wavey  jet  and  forms  the 
slipstream  typically  associated  with  propulsion  devices  operating 
in  a  fluid . 

The  previous  calculations  are  all  based  on  a  streamwise  ve¬ 
locity  U  at  the  foil  which  includes  the  induced  velocity  of  the 
slipstream  u  and  is  assumed  to  be  uniform.  In  order  to  determine 
the  free  stream  velocity  at  infinity  upstream  of  the  foil  it 
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is  necessary  to  compute  the  mean  value  of  u  or  u.  The  free  stream 
velocity  is  thus  taken  to  be 


U  =  U  -  u  [29] 

o 

The  instantaneous  induced  velocity  u  can  be  obtained  from  the 
induced  angle  a  caused  by  the  vortex  wake.  The  magnitude  of  the 
induced  velocity  is  taken  to  be  V  tan  a  and  is  directed  normal  to 
the  stream  direction  at  the  quarter  chord.  Thus  the  streamwise 
component  of  the  induced  velocity  is  given  by 


u  =  -  V  tan  sin  (a-Aa)  [24] 

The  induced  angle  0^  is  simply  the  difference  between  the 
instantaneous  kinematic  angle  of  attack  at  the  3/4  chord  and  the 
actual  angle  of  attack. 


a 


i 


a  ' 
e 


[25] 


The  last  term  is  the  steady  state  induced  angle  of  attack  due  to 
finite  span.  It  appears  here  because  the  Theodorsen  functions 
used  in  this  paper  have  been  normalized  with  respect  to  their 
steady  state  values.  If  the  Theodorsen  functions  were  not 
normalized,  then  the  two-dimensional  lift  curve  slope  would  have 
to  be  used  in  computing  the  forces  on  the  foil  and  all  three-di¬ 
mensional  effects  would  be  automatically  accounted  for  by  using 
the  three-dimensional  values  of  the  Theodorsen  functions. 
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The  mean  induced  velocity  can  now  be  expressed  as 


2mr 


u  =  -  7^;  j  V  tan  (a^)  sin  (a-Aa)d(cut) 


o 


[28] 


The  propulsive  efficiency  is,  of  course,  given  by 

TU 


n  = 


[31] 


where  T  and  P  are  obtained  from  Equations  [26]  and  [27] 


Computer  Calculations 

Numerical  calculations  of  the  preceding  theory  were  carried 
out  on  an  IBM  1130  digital  computer.  All  the  forces  and  moments 
acting  on  the  foil  were  computed  every  2.5  degrees  during  one 
cycle  of  foil  operation.  The  mean  values  were  obtained  by  a 
Simpsons  Rule  type  of  integration. 

The  values  of  the  function  C(k)  were  obtained  by  a  curve  fit 
to  the  coefficients  in  the  expression  obtained  by  Jones  for  el¬ 
liptic  wings  of  finite  span.  (Reference  4  )jones  obtained  an  ex¬ 
pression  of  the  form 


ik 


[  7A] 


C(k)  =  (F+iG)  =  1-Ci 


lk+C2 
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The  real  and  imaginary  parts  of  this  yield  expressions  for  P  and 

G. 


F  =  1  - 


Ci  kc 


U  =  - 


k2  +  C22 


CiC2k 


[?B] 


k2  +  C22 


J 


The  coefficients  Ci  and  C2  are  represented  by  the  expressions 


Ci 


C2 


0.5 


A 

A  +  2.32 


0.  l8l 


,  0.772 
A 


J 


[  7C] 


These  are  simply  empirical  curve  fits  to  the  theoretical  values 
of  Ci  and  C2 .  This  representation  gives  the  normalized  values 
of  F  and  G  when  k  is  based  on  the  average  chord  c  defined  as 

S 

C  5  b 

The  values  of  F  and  G  presented  in  Figures  2A  and  2B  are  obtained 
from  these  expressions. 

The  steady  state  foil  performance  was  represented  by  the 


following  expressions 


I 
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°L  ~  °L  “  H  " 

a 


[32] 


CL  =  CL  “  > 

max 


[33] 


C_,  =  CM  +  C..  a 
MM  M 
o  a 


[34] 


CD=CD  + 
o 


'CT  a 
1  La 


[35] 


where 


is  the  lift  coefficient. 


is  the  lift  curve  slope, 
a 

C  is  the  stall  lift  coefficient, 


is  the  1/4  chord  moment  coefficient, 
is  the  zero  lift  moment  coefficient. 


is  the  moment  curve  slope, 
is  the  drag  coefficient, 

is  the  zero  lift  drag  coefficient, 

is  the  geometric  aspect  ratio,  and 
is  the  efficiency  factor. 
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The  product  QE  is  interpreted  as  the  effective  aspect  ratio  in 
the  sense  that  the  induced  drag  is  the  same  as  for  an  elliptic 
wing  of  aspect  ratio  AE.  This  effective  aspect  ratio  was  used 
in  the  computation  of  the  function  C(k). 

The  previously  described  stall  criteria  were  imposed  to 
determine  when  unsteady  stall  would  occur.  If  stall  did  occur 
it  was  assumed  that  the  slipstream  velocity  would  be  reduced  by 
the  ratio  of  the  stalled  lift  coefficient  to  the  corresponding 
unstalled  lift  coefficient,  l.e.  C  /C  a. 

Xj  Ij 

a 

Numerical  calculations  were  carried  out  for  the  theoretical 
performance  of  the  model  tested  in  the  experimental  portion  of 
this  study.  These  are  presented  lauer  in  the  report  and  are  dis¬ 
cussed  in  conjunction  with  the  experimental  results. 

TEST  PROGRAM 

An  experimental  program  was  undertaken  to  determine  the  va¬ 
lid!  cy  of  the  previously  described  analysis.  The  tests  were  con¬ 
ducted  on  an  aspect  ratio  3  model  in  the  High  Speed  Channel  at 
HYDRONAUTICS,  Incorporated  and  cover  a  wide  range  of  operating 
conditi ons . 

Experimental  Apparatus 

Water  Channel  -  The  tests  were  conducted  in  the  High  Speed 
Water  Channel  at  HYDRONAUTICS,  Incorporated.  This  is  a  free 
surface,  circulating,  variable  pressure  water  channel  with  a 
test  section  2  feet  wide  and  12  feet  long.  A  sluice  gate  at  the 
entrance  of  the  test  section  allows  the  water  depth  to  be  varied 
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from  about  8  inches  to  2  feet.  The  present  tests  were  conducted 
at  atmospheric  pressure  with  a  water  depth  of  19  inches.  At  this 
depth  water  speeds  up  to  25  feet  per  second  can  be  obtained.  The 
maximum  Reynolds  number  obtained,,  based  on  chord  length  and  free 
stream  velocity,  was  thus  only  about  8  x  10  .  This  value  is 
rather  low  for  the  present  type  of  work  because  the  foil  stall 
characteristics  may  be  affected  at  this  Reynolds  number.  A  de¬ 
tailed  description  of  the  facility  used  is  presented  in  Reference  8 

Experimental  Model  -  A  model  with  a  4  inch  cnord  and  12  inch 
span  was  selected  as  a  good  compromise  between  a  model  large  enough 
to  obtain  reasonable  forces  and  yet  not  so  large  as  to  be  exces¬ 
sively  influenced  by  the  channel  boundaries.  The  foil  has  a  rec¬ 
tangular  planform  with  slightly  rounded  tips  and.  an  NACA  63AOI5 
section.  It  is  supported  horizontally  by  three  2  inch  chord  struts 
Two  of  these  struts  are  located  3  inches  from  the  foil  tips  and 
1-1/2  inches  from  the  foil  leading  edge  while  the  third  strut  is 
attached  to  a  pod  on  the  foil  centerline  and  located  4-1/2  inches 
from  the  leading  edge.  The  effective  foil  hinge  point  is  located 
midway  between  the  fore  and  aft  struts  at  the  3/^  chord.  The  for¬ 
ward  struts  were  located  inboard  of  the  foil  tips  so  that  they 
could  not  act  as  end  plates  and  thus  influence  the  sheding  of 
vorticity.  The  model,  at  mid-stroke,  was  located  10  inches  from 
the  channel  bottom  and  had  9  inches  of  water  above  it.  The  model 
is  shown  in  Figure  5. 
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Instrumentation  -  The  struts  supporting  the  foil  were  con¬ 
nected  through  variable  reluctance  force  gages  and  support  beams 
to  the  pistons  of  a  Planar  Motion  Mechanism  (PMM).  This  equipment 
is  shown  in  Figure  6.  These  pistons  of  the  PMM  are  in  turn  con¬ 
nected  to  cranks  so  that  the  forward  pair  of  struts  and  aft  strut 
can  be  oscillated  with  approximately  sinusoidal  motion.  The  system 
is  arranged  so  that  the  stroke  and  phase  relation  between  the  fore 
and  aft  struts  can  be  varied  through  a  wide  range.  The  resulting 
foil  motion*  referenced  to  the  J>/b  chord*  can  be  described  as  a 
combination  of  heave  h  and  angular  rotation  a  by  the  relations 

h  =  h  sin  (cat) 
o 

a  =  sin  (cot  +  0)  =  ai 

where  h  is  the  maximum  heave  displacement*  a  is  the  maximum 
angular  displacement*  and  0  is  the  phase  relation  between  the 
heave  and  angular  motion.  Since  the  mean  angle*  a^*  is  zero  in 
these  tests*  a  is  equal  to  ai. 

The  force  gages  are  arranged  to  measure  the  lift  and  drag 
forces  acting  on  the  forward  pair  of  struts  and  the  aft  strut.  In 
order  to  obtain  the  power  absorbed  by  the  foil  it  is  necessary  to 
multiply  the  instantaneous  lift  force  by  the  vertical  velocity  of 
the  corresponding  strut.  The  strut  velocities  were  obtained  from 
two  sin-cos  potentiometers  connected  to  the  PMM  drive  shaft  and 
set  in  phase  with  the  motion  of  each  piston.  Since  the  motion  of 
each  piston  is  essentially  sinusoidal  its  velocity  will  be  pro¬ 
portional  to  the  cosine  of  its  angular  displacement.  Thus  by 


HYDRONAUTICS ,  Incorporated 


-27- 

multiplying  the  lift  gage  signal  by  the  output  of  the  cosine  po¬ 
tentiometer  a  signal  proportional  to  power  was  obtained. 

Mean  values  of  thrust  and  power  were  obtained  by  feeding 
the  gage  signals  into  an  electrical  integration.  In  the  present 
tests  the  integration  was  carried  out  over  10  cycles.  This  ef¬ 
fectively  filters  out  the  hash  which  inevitably  appears  in  os¬ 
cillator  tests.  A  more  complete  description  of  the  Planar  Motion 
Mechanism  is  presented  in  Reference  8. 

Test  Procedure 

Steady  State  Performance  -  In  order  to  predict  the  oscillatory 
performance  from  the  theory  it  is  necessary  to  know  the  steady  state 
performance  of  the  foil.  This  could  be  predicted  from  theory,  how¬ 
ever  it  is  more  accurate  to  measure  this  performance  directly.  By 
making  these  measurements  with  the  foil  located  in  the  water  chan¬ 
nel  test  section,  the  influence  of  the  channel  boundaries  is  auto¬ 
matically  included  in  the  results.  It  has  been  assumed  that  the 
induced  velocities  from  the  image  of  the  vortex  system  in  the  chan¬ 
nel  boundaries  is  influenced  by  the  foj  1  oscillations  in  the  same 
manner  as  the  induced  velocities  from  the  primary  vortex  system. 

This  will  be  true  if  wavemaking  on  the  free  surface  can  be  ignored. 
The  present  tests  were  conducted  with  a  submergence  of  more  than 
two  chords.  Except  for  a  few  tests  at  zero  forward  speed,  the 
Proude  numbers  based  on  semichords  ranged  from  4.3  to  10.8.  Ac¬ 
cording.  to  the  results  of  Reference  1,  there  should  be  no  appre¬ 
ciable  influence  of  the  free  surface  under  these  conditions  for 
the  range  of  reduced  frequencies  involved.  Froude  numbers  based 
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on  water  depth  ranged  from  1.4  to  3*5  with  the  exception  of  the 
zero  speed  tests  where  the  Froude  number  is,  of  course,  zero.  The 
channel  flow  is  thus  sufficiently  super  critical  that  wavemaking 
effects  due  to  finite  depth  will  be  negligible  except  for  the  zero 
speed  tests.  Therefore.,  no  additional  boundary  corrections  should 
be  required  in  predicting  the  performance  of  the  oscillating  foil 
in  the  water  channel  when  the  predictions  are  based  on  the  steady 
state  foil  performance  as  measured  in  the  water  channel. 

Results  of  the  steady  state  performance  measurements  are  pre¬ 
sented  in  Figure  7.  They  were  obtained  at  speeds  of  13  and  23  feet 
per  second  and  cover  an  angle  of  attack  range  from  -26  to  +26  de¬ 
grees.  The  coefficients  obtained  from  these  tests,  and  used  in 
the  theoretical  predictions,  are  listed  below.  Pitching  moment 
measurements  were  not  made.  Because  the  foil  is  symmetric  it  was 
assumed  that  the  quarter  chord  pitching  moment  would  be  zero. 

C  =  0.065  per  degree 
a 

CL  =  0.90 

max 

cD  =  0.027 

o 

E  =1.05 

The  influence  of  the  channel  boundaries  is  reflected  in  these 
results.  The  measured  lift  curve  slope  Is  about  20  percent  higher 
than  would  be  expected  in  open  water  and  the  value  of  E  being 
greater  than  one  Indicates  that  the  induced  drag  is  the  same  as  an 
elliptic  wing  of  slightly  higher  aspect  ratio.  We  would  expect 
E  to  have  a  value  of  about  0.95  in  open  water.  The  performance 
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used  in  the  predictions.,  based  on  these  coefficients,  is  Indi¬ 
cated  by  the  dashed  lines  in  Figure  7- 

Oscillating  Tests  -  Before  proceeding  with  a  discussion  of 
the  oscillating  tests  it  Is  necessary  to  define  the  parameters 
and  nondimensional  coefficients  used  in  presenting  the  results. 

As  previously  mentioned  the  foil  motion,  referenced  to  the  3/^  chord 
can  be  described  as  a  combination  of  heave  h  and  angular  rotation  a 
by  the  relations 

h  =  n  sin  (27rnt) 
o 

a  =  a  sin  (27rnt  +  Q) 
o 

where  the  oscillating  frequency  od  in  radians  per  second  has  been 
replaced  by  the  frequency  in  cycles  per  second  n.  The  independent 
test  parameters  are  listed  below 

n  oscillating  frequency, 

h  maximum  heave  displacement, 
o 

aQ  maximum  angular  displacement, 

6  phase  relation  between  heave  and  angular  displace¬ 
ment,  and 

Uq  free  stream  velocity. 

The  quantities  measured  are,  of  course, 

T  the  mean  thrust  delivered,  and 

p  the  mean  power  absorbed. 


[36B] 


,  J 
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Two  sets  of  nondimensional  coefficients  can  be  defined.  The 
first  of  these  sets  is  analogous  to  the  usual  thrust,  power,  ad¬ 
vance  coefficient,  and  blade  area  ratio  used  in  presenting  con¬ 
ventional  propeller  data  and  are  defined  below  for  an  oscillating 
foil. 


KT 


_ T _ 

pn2(2ho)2(2ho  x  b) 


KP 


_ P _ 

pn3(2ho)3(2ho  x  b) 


J 


u 


2h  n 
o 


BAR 


S 

2h  x  b 
o 


These  coefficients  are  useful  in  comparing  the  performance  of  an 
oscillating  foil  propulsor  with  a  conventional  propeller  but  be¬ 
come  undefined  for  the  case  of  h  =0. 

o 

A  second  set  of  coefficients  more  suited  to  the  oscillating 
foil  propulsor  can  be  defined  as: 


1 
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"p 

CKP  =  — £ - 

pn3c3S 

U 

j.  =  -2.  =  2L 

nc  k 

h  *  =  h  /c 
o  o/ 

These  coefficients  are  used  in  the  presentation  of  the  experi¬ 
mental  data  in  this  report. 

The  propulsive  efficiency  rj  is  defined  in  the  usual  manner 
as 

T-  U 

o  _  KT-J  CKT-J'  r  n 

R  KP  CKP  l3U 

r 

The  relation  between  the  two  sets  of  coefficients  is  given  below 

KT  =  CKT-(BAR)3 
KP  =  CKP- (BAR)4 
J  =  J'/2h 

Optimum  performance  of  an  oscillating  foil  propulsion  sys¬ 
tem  will  be  obtained  when  the  maximum  transverse  velocity  of  the 

foil  is  of  the  same  order  as  the  free  stream  velocity.  This  cor- 

* 

responds  to  J  =  t  or  J 1  =  2-rrh^  .  In  the  present  tests  physical 
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limitations  of  the  test  equipment  restricted  the  maximum  attain¬ 
able  transverse  velocity  to  0.628  of  the  free  stream  velocity 
(J  =  5»0)>  This  was  obtained  with  the  maximum  oscillating  fre¬ 
quency  of  5  cps,  the  maximum  stroke  of  0.2  feet  and  the  minimum 
water  speed  at  which  wavemaking  could  be  neglected  of  10  fps. 

j#. 

These  correspond  to  h^  =0.6  and  J'  =  6. 

It  was  desired  to  obtain  data  on  the  influence  of  as  many 
of  the  various  parameters  as  possible.  The  range  of  test  param¬ 
eters  covered  is  shown  in  Table  I.  The  experimental  and  theo¬ 
retical  values  of  CKT,  CKP,  and  r\  have  been  plotted  as  functions 
of  a ,  Q,  and  J'.  Tables  II,  III,  and  IV  provide  a  summary  of  the 
data  obtained  and  list  the  figure  numbers  where  each  parameter  is 
presented.  In  addition  to  the  performance  data,  these  figures 
include  the  theoretical  estimate  of  the  fraction  of  each  cycle 
that  the  foil  is  stalled. 

Accuracy  -  Because  of  the  oscillatory  nature  of  the  measured 
forces,  the  maximum  loads  experienced  by  the  force  gages  is  con¬ 
siderably  larger  than  the  mean  values  recorded.  This  effect  is 
particularly  important  in  measurements  of  the  mean  power  because 
the  inertial  forces  on  the  model  do  not  contribute  to  the  mean 
»  value  of  the  absorbed  power.  The  force  gages  must,  of  course,  be 
capable  of  measuring  the  peak  loads.  As  a  result,  the  accuracy 
of  the  dynamic  measurements  is  somewhat  less  than  that  obtainable 
in  static  tests.  The  accuracy  of  the  measured  quantities  and  the 
corresponding  coefficients  is  estimated  to  be  within  the  following 
limits : 


t 


HYDRONAUTICS ,  Incorporated 


-33- 

Lift  (200  lbs  maximum) 
Drag  (25  lbs  maximum) 
Dynamic  pressure 

(4.4  psi  maximum ) 
Velocity  (25  fps  maximum) 
Oscillating  frequency 
(5  cps  maximum) 
a 

o 

0 

h  (0.2  feet  maximum) 

CKT  (  5  maximum) 

CKP  (260  maximum) 

J '  (19  maximum ) 

h  *  (0.6  maximum) 
o  ' 


±2.0  lb 
±0.25  lb 

±0.1  psi 
±0.03  fps 

± 0.  05  cps 
±0.5° 

±1.0° 

±0. 005  feet 

±0.  2 
±4.  2 
±0. 05 
±0.  015 


Results  and  Discussion 

An  example  of  the  predicted  forces  and  moments  on  the 
model  during  one  cycle  of  oscillation  is  presented  in  Fig¬ 
ure  8.  The  particular  combination  of  parameters  in  this  ex¬ 
ample  yields  a  relatively  high  efficiency  (57  percent)  for 
the  thrust  coefficient  and  speed  involved.  The  foil  is 
stalled  44  percent  of  the  cycle  which  results  in  the  slightly 
flattened  portion  of  the  lateral  force  curve  (CKS)  and  moment 
curve  (CKQ).  The  influence  of  the  stall  can  also  be  seen  in 
the  thrust  curve  (CKT).  If  stall  did  not  occur,  the  thrust 
curve  would  have  a  nearly  sinusoidal  shape.  In  the  percent 
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case,  angular  rotation  of  the  foil  caused  stall  to  occur  at  a 
lift  coefficient  of  0.5  while  under  steady  conditions  the  foil 
would  reach  a  lift  coefficient  of  0.9  before  it  stalled. 

Figure  9  shows  the  predicted  performance  of  the  model  foil 

— 

as  a  function  of  a  and  9  for  the  case  of  h  =0.6  and  J'  =9* 

o  o 

This  figure  shows  that  the  maximum  thrust  occurs  at  a  phase  angle 
of  about  60  degrees  while  the  maximum  efficiency  occurs  at  a 
phase  angle  of  about  75  degrees.  At  higher  values  of  J'  the  maxi¬ 
mum  thrust  and  efficiency  occur  at  lower  phase  angles  while  at 
low  values  of  J1  the  reverse  is  true.  For  example,  at  a  J'  of 
15  the  maximum  thrust  and  efficiency  occur  at  phase  angles  of 
about  20  and  40  degrees  respectively. 

Figures  10  through  66  show  both  the  theoretical  and  the  ex¬ 
perimental  performance  of  the  foil.  These  figures  indicate  that, 
with  the  exception  of  the  zero  forward  speed  data  (Figure  28)  the 
theory  is  in  general  agreement  with  the  data.  However,  there  are 
several  trends  which  can  be  noted.  At  high  values  of  J'  and  low 
values  of  9  the  theory  tends  to  overestimate  CKT  while  at  high 
values  of  6  the  theory  overestimates  CKT  very  slightly.  At  low 
values  of  J',  particularly  when  6  is  also  small,  the  theory  under¬ 
estimates  both  CKT  and  r\.  There  is  also  a  general  trend  for  the 
peak  thrust  and  peak  efficiency  to  occ,  at  slightly  lower  values 
of  9  than  predicted. 

At  zero  advance  coefficient  the  theory  predicts  very  low 
levels  of  thrust  for  the  oscillating  foil  while  this  is  the  con¬ 
dition  where  a  conventional  propeller  usually  reaches  its  maximum 
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thrust.  The  reason  for  this  can  be  seen  by  comparing  the  advance 
coefficient  for  the  two  systems  at  the  normal  operating  condition. 

A  conventional  propeller  normally  operates  with  an  average  blade 
velocity  in  a  direction  normal  to  the  free  stream  that  is  con¬ 
siderably  in  excess  of  the  free  stream  velocity  (J  <  ir) .  When 
this  propeller  is  brought  to  a  condition  of  zero  advance  the  dy¬ 
namic  pressure  is  only  slightly  reduced  while  the  blade  angle  of 
attack  is  increased.  The  increased  lift  coefficient  resulting 
from  the  increased  angle  of  attack  is  normally  more  than  sufficient 
to  offset  the  reduction  in  dynamic  pressure.  The  result  is  an  in¬ 
crease  in  thrust. 

The  present  oscillating  foil  with  hQ*  =0.6  operates  ef¬ 
ficiently  with  transverse  velocities  considerably  below  the  free 
stream  velocity  (J  >  tt).  Thus  when  the  foil  is  brought  to  zero 
advance  nearly  all  the  dynamic  pressure  is  lost.  The  resulting 
angles  of  attack  are  so  large  that  the  foil  is  stalled  nearly  the 
entire  cycle.  The  result  is  a  large  loss  In  thrust. 

The  theory  indicates  that,  under  these  conditions,  stall 
cannot  be  avoided  by  increasing  to  reduce  the  angle  of  attack. 
Increasing  increases  the  angular  velocity  (increased  y)  which 
also  causes  stall.  Stall  can,  however,  be  reduced  at  zero  ad¬ 
vance  by  using  a  sufficiently  large  value  of  h  *.  If  adequately 

-fc- 

large  values  of  h^  are  selected  the  thrust  can  be  made  to  in¬ 
crease  with  decreasing  advance  ratios  in  a  manner  similar  to  a 
conventional  propeller.  In  general,  the  theory  Indicates  best 
performance  will  be  obtained  with  the  largest  practical  values 

of  h  *. 
o 
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It  can  be  seen  in  Figure  28  that  the  model  foil  produces  con¬ 
siderably  higher  thrust  at  zero  advance  than  the  theory  predicted. 
The  assumptions  underlying  the  analysis  are,  of  course,  not  met 
under  these  conditions.  The  test  results  rnay  also  have  been  in¬ 
fluenced  by  wavemaking  which  is  not  accounted  for  in  the  analysis. 
The  actual  cause  of  the  discrepancy  is  not  adequately  understood. 

It  should  be  pointed  out  that  the  maximum  potential  efficiency 
of  an  oscl  1 ating  foil  propulsor  is  not  as  high  as  for  a  conven¬ 
tional  rotating  propellc.  of  the  same  "disc"  loading.  The  fluc¬ 
tuating  loading  on  the  oscillating  foil  results  in  shedding  of 
vorticity  in  the  wake  beyond  that  of  the  trailing  vortex  system 
associated  with  a  steady,,  finite-span  wing.  This  represents  an 
additional  loss  of  energy  which  is  accounted  for  by  the  Theodorsen 
functions.  The  disc  loading  (thrust  loading  over  the  area  swept 
by  the  foil)  is  also  not  uniform.  This  results  in  additional  mo¬ 
mentum  losses  associated  with  a  nonuniform  slipstream.  Finally, 
because  the  foil  must  operate  over  a  wide  range  of  lift  coeffi¬ 
cients,  it  can  not  be  "optimiz-d"  for  a  specific  lift  coefficient. 
This  results  in  both  a  reduction  in  the  average  lift-drag  ratio 
daring  one  cycle  and  a  reduction  in  the  mean  thrust  that  can  be 
produced  without  cavitation  or  blade  stall.  However,  it  will 
often  be  possible  to  obtain  a  greater  swept  area  with  the  oscil¬ 
lating  foil  than  with  the  conventional  prcpell-r.  The  increased 
efficiency  thus  obtained  may  well  be  sufficient  to  give  the  os¬ 
cillating  foil  propulsor  higher  efficiency.. 
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FOIL  SYSTEM  FOR  "SKI  BARGE" 

The  previously  described  theoretical  analyses  were  used  for 
the  preliminary  design  of  an  oscillating  foil  propulsor  mutable 
for  use  on  a  15  knot,  shallow-draft  (eight  inches)  boat  of  2000 
pound  payload,  such  as  the  "Ski  Barge."  A  variety  of  foil  geo¬ 
metries  and  motion  parameters  were  examined  to  determine  a  suitable 
compromise  between  requirements  for  high  static  thrust,  high  cruise 
efficiency,  and  minimum  blade  stalling  over  the  operating  range. 
Blade  stall  is  of  particular  importance  because  the  separated  flow 
may  lead  to  problems  of  foil  ventilation  for  a  system  of  such 
shallow  immersion. 

The  foil  arrangement  selected  is  shown  in  Figure  r.7,  It  con¬ 
sists  of  two  groups  of  three  foils  each.  Each  group  is  mounted 
on  the  after  edge  of  a  parallelogram  frame,  the  trailing  arms  of 
which  are  15  inches  long.  These  arms  are  swung  laterally  through 
±45  degrees.  The  two  groups  of  foils  are  oscillated  in  opposite 
directions  so  that  no  lateral  vibration  is  transmitted  to  the 
stern  of  the  boat.  The  foils  are  8  inches  deep  and  have  a  mean 
chord  of  6  inches.  They  are  mounted  between  cheek  plates  in  such 
a  way  that  they  can  fold  aft  if  the  boat  is  grounded.  The  cheek 
plates  hinge  about  a  vertical  axis  that  is  3  inches  forward  of 
the  mean  quarter  chord  and  have  a  ventilation  fence  at  their 
lower  end.  The  foil  was  assumed  to  have  a  lift  curve  slope 

(C  )  of  0.033  Per  degree,  a  zero-lift,  drag  coefficient  (C  ) 
a  o 

of  0.012  ana  a  stall,  lift  coefficient  of  0.90. 
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Best  all  around  performance  was  obtained  with  a  value  of 
aQ  of  40  degrees,  0  of  90  degrees,  and  an  oscillating  frequency 
of  360  RPM.  The  calculated  performance  of  one  foil  is  shown  in 
Figure  68  for  speeds  up  to  20  knots.  Relatively  high  efficiency 
was  obtained  in  spite  of  the  low  aspect  ratio  foils  because  the 
large  swept  area  produces  rather  low  "disc"  loading.  At  zero  for¬ 
ward  speed  the  system  will  produce  1000  pounds  of  thrust  with 
50  horsepower.  Under  these  conditions  the  foil  would  be  stalled 
74  percent  of  each  cycle.  Foil  stall  would  stop  at  8  knots  and 
at  15  knots  the  system  would  deliver  570  lbs  of  thrust  with  4l 
horsepower  at  an  efficiency  of  64  percent.  Slightly  better  per¬ 
formance  could  be  obtained  at  both  high  and  low  speeds,  if  con¬ 
trol  were  provided  over  and  0.  This  would  be  similar  to  the 
performance  gains  obtainable  with  a  controllable  pitch  propeller. 
Of  course,  higher  thrust  levels  can  be  obtained  with  higher  in¬ 
stalled  power  but  the  efficiency  will  be  rapidly  reduced  and  in¬ 
creased  stalling  will  occur  at  the  lower  speeds. 

Steering  can  be  accomplished  by  supplying  a  mean  angle  of 
attack  to  the  foils  in  a  manner  similar  to  a  rudder.  This  will 
produce  a  lateral  force  vector  even  at  zero  forward  speed.  In 
fact,  if  complete  control  over  the  motion  parameters  is  supplied, 
this  type  of  foil  system  will  produce  thrust  vectors  in  any  direc¬ 
tion  as  well  as  yawing  moments. 

The  foregoing  analyses  show  that  an  oscillating  foil  pro- 
pulsor  can  provide  efficient  shallow  water  propulsion  with  a  high 
degree  of  maneuverability.  Its  ultimate  practicality  will  depend 
more  on  the  solution  of  mechanical  drive  problems  than  on  the 
hydrodynamic  performance. 
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|  Y5PV2Sc 


[5] 


[  (-r-  +  d)  a  sin  a  +  U  +  i]  +  [  (  -r-  +  d )  d  cos  a-ft] 


C(k)  =  P(k)  +  iG(k)  =  1-Ci 


[6] 

[7A] 


F  =  1  - 


Ci  k 


2  *7 


k2+  C22 


G  =  - 


CiC2k 

k2+  C22  ^ 


[7B] 


Ci  =  0.5 


A 


A+2. 32 


*7 


c2  =  0.181  +  °'J7-  j 


[7C] 


k  = 

00  c/2 

U 

[8] 

a  = 

a  +  ai 
m 

[9] 

a 

=  a  +  a  1 

[10] 

e 

m  e 

a  ' 
e 

'  S=A  Vp2  +  GS 

[11] 
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a^c/n  -  a i  +  arctan 

(c/g  + 

d)  a  cos  ct  -  ft 

[12] 

_(c/2  +  d)  a  sin  a  +  'J+  £ 

a3c/4  =  a3c/4  evaluated  at  time  (t-At) 

[13] 

arctan  (-G/F) 

^  O) 

[14] 

Act  =  -  27 

[15] 

—  — 

,  c  a 

7  =  arctan  — 

[16] 

wc  =  (U  +  £ )  cos  a  + 

h  sin  a 

[17] 

Mass  Forces 

LQm  =  L2  cos  a  -  D2  sin  a 

DQm  =  D2  cos  a  +  L2  sin  a 

> 

[18] 

Mam  =  M2  -  ( c/4  +  d )L2 

L2  =  am*w_T 

N 

D2  =  am • (t/c  )2  w c 

[19] 

M2  =  -am *K(c/4 )2-  a 
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am  = 


Vl+A2 


P  ^ ( c/2 ) 2  b 


[20] 


=  [2  +  h  &]  sin  a  +  [(U+  i  )a-h  ]  cos  a  +  [c/4  +  d]a 


wp  =  [ £+  ha]  cos  a  -  [ ( U +  H  )a  -  h]  sin  a 


[21] 


r 


K  = 


0.30 

{  0. 40 

0.33 


A  — ►  00 
A  =  1 
A  0 


[  22A  ] 


K  ~  A+1.43  +  °-33 


[  22B] 


Foil  Mass  Inertial  Forces 


L  =  m  [(c/4  +  d-s)(a  cos  a -a2 sin  a)  -  h] 


D  =  m  [(c/4  +  d-s ) (a  sin  a  +  a2 cos  a)  +  iH 
m  ' 


M^  =  171  ^  +  4l-s)(h  cos  a  -  l  sin  a)- 


)  [23] 


{V2  +  (c/4  +  d-s  )2  a  ]  J 
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u  =  -  V  tan  a.  sin  (a-Aa) 
1  '  ' 


a . 
1 


“3c/4  -  v  +  °lAa 


2nrr 


T  =  - 


2mr 


D  d(ajt) 


o 


2ri7r 


P  =  - 


2ntr 


(Lh-D i+  Met)  d(aot) 


2n7r 


u  = 


2mr 


V  tan  (a^  )sin  (a-Aa  )d  (ait ) 


U  =  U  -  u 
o 


a)  c/2 


HYDRONAUTI CS ,  Incorporated 


-44- 

APPENDIX  II 

DERIVATION  OP  EQUATIONS  FOR  THE 
HYDRODYNAMIC  FORCES  AND  MOMENTS 

The  equations  for  the  hydrodynamic  forces  and  moments  acting 
on  an  oscillating  foil  have  been  derived  from  those  developed  by 
Theodorsen  (Reference  6)  and  Garrick  (Reference  3).  Theodorsen 
treated  the  problem  of  a  thin,,  uncambered  foil  and  obtained  a 
linearized  solution  for  the  case  where  the  oscillatory  velocities 
are  small  compared  to  the  free  stream  velocity.  It  is  important 
to  note  that  the  lift  L  and  moment  M  are  of  order  a  (first  order) 
while  the  drag  D  (or  thrust)  is  of  order  La  or  a  squared  (second 
order).  Because  the  drag  terms  are  of  higher  order  they  were 
dropped  in  Theodorsen' s  linear  theory.  The  wake  thus  computed 
from  this  linearized  theory  contains  only  the  induced  velocities 
which  are  normal  to  the  free  stream  direction  and  does  not  con¬ 
tain  Induced  velocities  in  the  streamwise  direction. 

Garrick  has  computed  second  order  terms  for  the  forces  on 
the  foil  arising  from  the  so-called  leading  edge  suction.  This 
force  is  proportional  to  (S')2  where  S'  is  given  in  Equation 
4.1-AII  below.  However,,  he  did  not  compute  the  induced  veloc¬ 
ities  in  the  streamwise  direction.  As  a  result  there  is  no 
slipstream  in  his  analysis.  This  yields  an  optimistic  propul¬ 
sive  efficiency  since  the  kinetic  energy  lost  in  the  slipstream 
has  not  been  accounted  for. 

If  we  let  the  induced  slipstream  velocity  at  the  foil  be  de¬ 
signated  u  and  the  free  stream  velocity  at  infinity  upstream  of 
the  foil  be  U  .  Then  the  total  streamwise  velocity  U  at  the 
foil  is 
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U  =  U  +  u  [ 1.0A-AII] 

o 

Since  it  is  not  possible  to  compute  the  slipstream  velocity 
at  the  foil  until  the  foil  loading  is  known^  and  visa-versa,  one 
of  these  two  quantities  must  be  assumed.  The  approach  taken  here 
is  to  compute  the  forces  and  moments  on  the  foil  based  on  the 
total  local  velocity  U  and  then  compute  the  induced  velocity  u. 

The  mean  value  of  u  is  then  computed  and  the  free  stream  velocity 
is  taken  to  be 

Uq  -  U  -  u  [1.0B-AII] 

where  the  bar  indicates  average  values. 

The  following  equations  for  lift,  drag,  and  moment  on  a 
thin,  uncambered,  oscillating  foil  are  taken  from  Reference  3  and 
have  been  converted  to  the  notation  of  the  present  paper.  The 
reduced  frequency  k  is  based  on  the  local  velocity  U  in  the  stream- 
wise  direction  • 

L  =  pU27r(c/2)  [C(k)Q]  +  pir  ( c/2  )2  [  Ua-h  +  (c/4  +  d)a] 

[ 1 . 1 -All ] 


D  =  La  -  pm ( c/2 ) (S 1 )2 


[  2. 1-AII] 
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M  =  -  pU27r(c/2)d[  C(k)Q] 


p7r(  c/2  )2  [  ( c/2  +  d  )Ua  + 


(c/4)2 


+  (c/4  +  d)2 


a  -  (c/4  +  d  )h  ] 


[3.1-AIl] 


where 


S  ' 


[  2 c(k )Q  -  c/2  d] 


[ 4. 1-AII] 


Q  =  Ua  -  h  +  (c/2  +  d )  a 


[5. 1-AII] 


Equation  5. 1-AII  can  be  rewritten  in  terms  of  an  angle  by  dividing 
through  by  U 


Q 

U 


a 


h  (c/2  +  d)  a 
U  +  U 


[5.2-All] 


Note  that  this  is  the  kinematic  angle  of  attack  at  the  3/4  chord 
and  that  this  quantity  always  occurs  in  conjunction  with  the 
Theodorsen  function  C(k).  Therefore  we  will  define  an  "equivalent 
steady-state  angle  of  attack"  ag  as 

ag  =  C(k)  § 


Substituting  and  S'  into  the  equations  for  lift,  drag,  and 
moment  we  obtain 
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L  =  |pU2c  2ma  +  pm(  c/2  )2  [  Uct  -  h+(c/4  +  d)  a]  [1.2-AIl] 
D  =  La  -  |-pU2 c  2tt  ^[a  -  c/2  ~]2  [2.2-All] 

M  =  -  ^pU2c  2-ma^d  -  p7r(c/2)2[  (c/2  +  d  )U  a 

+  j^i(c/4)2  4-  (c/4  +  d)2  J  a  -  (c/4  +  d)  h] 

[3.2-AIl] 

The  equations  for  lift,  drag,  and  moment  are  now  expanded 
and  their  terms  grouped  so  that  they  can  be  identified 

L  =  [ipU2]  c  [2ir]ae  +  [  pm ( c/2  )2  ]  [  Ua  -  h  +  (c/4  +  d)  a] 

L  -  [|pU2]  S[CL  ]  ag  -  [am]  [wN] 

a 

L  =  Lx  +  L2  [1.3-AII] 

f>  =  La  -  [ipU2]  c  [2ir]ae[ae  -  2 (c/4  ^)] 

-  [|pU2]  c  [2m]  (  f  |)2 
D  =  (Li  +  L2 )  a  -  Li  (a  -  27 )  -  Si 


D  =  Li  (a-Aa)  -  Si  +  L2CI 


[  2. 3-AII] 
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M  =  -  [ipU2]c[27r]aed  -  [tpU2]c2[7r/2]  (  |  ?j) 

-  [p7r(c/2)2][  (c/4  +  d)(Ua  -  h  +  (c/4  +  d)  a) 
+  *(c/4)  a] 


M  =  -  Lid  -  L2(c/4  +  d)  +  -  am-K(c/4)2  a]  [3.3-All] 

where  the  following  quantities  are  identified 


=  2t r  is  the  two-dimensional  lift-curve  slope, 
a 

S  =  c  x  1  is  the  foil  area  per  unit  span, 

am  =  p7r(c/2)2  is  the  theoretical  two-dimensional  added  mass. 


a 

e 


[Ua  -  h  +  (c/4  +  d)  a]  is  the  acceleration  of  the  foil 

midchord,  normal  to  the  plane  of 
the  chord,  relative  to  the  fluid. 


C(k) 


a 


4  .  ( c/2  +  d )  a 

U  U 


is  the  instantaneous  angle 
of  attack  at  the  3/4  chord 
and  includes  the  unsteady 
effects  of  the  wake.  The 
effect  of  the  wake  is  to 
reduce  the  angle  of  attack 


by  the  factor  ~\/  F2  +  G2 
and  to  cause  a  time  lag  of 
At  =  [arctan  (-G/F)]/cd. 
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That  is 

a  =  a  -  ^  +  (c/2  +  ^  “\ / F2  +  G2  evaluated  at  time 

e  U  U  V 

J  ( t-At ) , 

Aa  =  ag  -  2y  is  the  instantaneous  angle  of  attack  at  the 
1/4  chord, 

c  ct 

y  =  —  is  the  change  in  angle  of  attack  per  1/4  chord. 

This  is  a  measure  of  the  effective  camber  generated 
by  the  foil  angular  velocity, 

The  following  forces  can  now  be  defined 

Li  =  fpU2S  C  a  is  the  instantaneous  foil  lift  due  to 
a 

hydrodynamic  circulation.  It  is  equal 

to  the  steady-state  lift  evaluated  at 

a  =  a  . 
e 

L2  =  am.w^  is  the  instantaneous  foil  lift  due  to  added 
mass. 

Si  =  ipU2S  2iry2  is  a  thrust  force  attributed  to  "leading 

edge  suction"  which  results  from  the 
foil  rotation, 

=  -|pU2S  c  —  y  is  a  pitching  moment  due  to  foil  rota¬ 
tion.  This  moment  is  the  same  as  that 
of  a  circular  arc  camber  foil  of  camber 
chord  ratio  equal  to  2y , 
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M2  =  -  am-K(c/4)2  a  is  the  pitching  moment  due  to  added  mass. 

The  product  [am  K(c/4)2]  can  be  recognized  as  the  moment  of 
inertia  of  the  added  mass  about  the  foil  midchord.  The  effective 
radius  of  gyration  of  the  added  mass  is  therefore  VI  (c/4)  where 
K  is  equal  to  1/2  for  a  two-dimensional  foil. 

The  preceding  equations  are  those  derived  from  small  amplitude 
potential  flow  theory  for  a  thin.,  uncambered,  foil  oscillating 
about  a  mean  angle  of  attack  of  zero.  In  order  to  include  viscous 
drag,  and  permit  cambered  foils  of  finite  thickness  oscillating 
about  a  mean  angle  of  attack  other  than  zero,  the  following  quan¬ 
tities  are  introduced  and  defined. 

Di  =  -|pU2S  CQ  is  the  foil  drag  associated  with  Li . 

is  the  steady-state  drag  coefficient 

evaluated  at  a  =  a  , 

e 

Mi  =  -|pU2Sc  CM  is  the  foil  quarter-chord  pitching 

moment  associated  with  Li .  C  is 

c/4 

the  steady-state  quarter-chord  pitching 

moment  coefficient  evaluated  at  a  =  a  , 

e 

D2  =  am(t/c)2  wc  is  the  instantaneous  foil  drag  due  to 

added  mass.  Here  the  effective  added 

mass  acting  in  the  chordwise  direction 

is  assumed  to  be  (t/c)2  times  the 

added  mass  acting  normal  to  the  chord. 

w  is  the  acceleration  of  the  foil  mid- 
0 

chord  in  the  plane  of  the  chord,  rela¬ 


tive  to  the  fluid. 
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is  the  mean  angle  about  which  the  foil  oscillates. 

The  Theodorsen  function  does  not  apply  to  a  since 

o 

the  oscillating  wake  does  not  effect  the  mean  angle 
of  attack.  The  equivalent  steady-state  angle  of 
attack  then  becomes 


a  =  a  +  C(k)  a  -  ^ 
e  o  U  U 


The  equations  for  lift,  drag,  and  moment  now  become 


L  =  Li  +  L12 


[ 1 . 4 -All ] 


D  =  Di  +  Li  (a-Aa)  -  S; 


[2.4-AIl] 


M  =  Mi  -  Lid  +  +  M2 


[ 3 • 4-AII ] 


The  mean  thrust  T  and  mean  power  P  can  be  determined  by 
integrating  the  forces  over  an  integer  number  of  cycles 


T  =  D  d(cot) 

2tt 


[4.4-AII] 


P  =  —  /  (L  ft  +  Ma)d(oyt) 

2t t  I 


[ 5 . 4-AII ] 
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The  large  amplitude  approximations  are  made  by  replacing 
linearized  values  of  the  various  quantities  by  their  actual  values. 
The  small  amplitude  theory  assumes  that  the  angles  are  small  and 
that  perturbation  velocities  are  small  compared  to  the  free  stream 
velocity.  Thus  the  sine  and  tangent  of  angles  are  taken  equal  to 
the  angle,,  the  cosine  is  taken  as  unity,  and  higher  order  perturba¬ 
tion  velocity  terms  are  neglected.  Many  of  the  important  non- 
linearities  can  be  recovered  by  simply  replacing  the  linearized 
velocity  terms  with  those  obtained  from  a  kinematic  analysis  of 
the  motion  and  by  including  the  correct  trigonometric  relations  in 
the  resolution  of  force  and  velocity  vectors.  The  resulting  equa¬ 
tions  are  those  which  appear  in  Appendix  I. 

While  these  equations  were  derived  for  a  foil  of  infiniet 
span  they  are  also  valid  for  foils  of  finite  span.  All  that  is 
required  is  that  the  two-dimensional  values  of  certain  coefficients 
be  replaced  by  their  corresponding  three-dimensional  values.  These 

include  C  ,  C  .  ,  C(k),  Ki  and  am.  Note  that  C  will  include 

^  u  c/4  U 

both  the  paracite  drag  and  the  induced  drag  that  would  normally 

occur  in  steady-state  performance  when  a  =  ag.  This  latter  re¬ 
sult  is  shown  to  be  correct  in  Reference  4. 

The  induced  velocity  u  can  be  obtained  from  the  induced 
angle  a^.  The  magnitude  of  the  induced  velocity  is  taken  to  be 
Uo^  and  is  directed  normal  to  the  stream  direction  at  the  quarter 
chord.  Thus  the  streamwise  component  of  the  induced  velocity  is 
given  by 


u  =  -  Uai(a-Aa) 


[  6-AII] 
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The  induced  angle  a,  is,of  coi'rsej  due  to  the  vortex  wake 
and  is  simply  the  difference  between  the  instantaneous  kinematic 
angle  of  attack  at  the  3/^  chord  and  the  actual  angle  of  attack 


a . 
1 


h  ,  (c/2  +  d)a 

u  +  u 


a. 

i 


V1^ 


[7-AII] 


The  last  term  is  the  steady  state  induced  angle  of  attack  due  to 
finite  span.  It  appears  here  because  the  Theodorsen  functions 
used  in  this  paper  have  been  normalized  with  respect  to  their 
steady  state  values.  If  the  Theodorsen  functions  were  not  nor¬ 
malized,  then  the  two-dimensional  lift  curve  slope  would  have  to 
be  used  in  computing  the  forces  on  the  foil  and  all  three-dimen¬ 
sional  effects  would  be  automatically  accounted  for  by  using  the 
three-dimensional  values  of  the  Theodorsen  functions. 


The  mean  induced  velocity  has  been  calculated  for  the  case 
of  a  two-dimensional  foil  in  pure  heave  (aQ  =  0)  and  can  be  ex¬ 
pressed  as 


2tt 


u 

U 


l-C(k) 


h 

U 


C(k) 


cot 


[8.1-AII] 


-  =  2k2 
U 


(h  *)2 
o 


F-(F2  +  G2) 


[8.2-AII] 
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where 


h  =  h  a)  sin  (cut ) ,  and 


h  *  =  h  /c. 
o  o7 


Thus 


U 

o 


U 


r  i-2ka(ho*)2 


P-  (F2 


[9-AII] 


For  the  case  of  pure  heave  the  thrust  is  positive  for  all  posi¬ 
tive  k  and  zero  when  k  =  0.  Prom  the  known  behavior  of  the 
Theodorsen  functions,  it  can  be  shown  that  P  >  (F2  +  G2 )  for 
k  >  0  and  the  equality  holds  only  if  k  =  0.  Therefore  the  induced 
velocity  is  always  positive  under  conditions  of  positive  thrust 
and  zero  when  the  thrust  is  zero,  as  expected.  The  free  stream 
velocity  at  infinity  upstream  of  the  foil  is,  of  course,  always 
less  than  the  velocity  at  the  foil.  The  ratio  of  free  stream 
velocity  to  local  velocity  (U  /U)  is  shown  in  Figure  69  for  values 

of  h  *  of  0,  0.1,  1.0,  and  10. 
o 

It  can  be  seen  from  the  preceding  equations  that  the  slip¬ 
stream  is  a  second  order  function  of  both  the  reduced  frequency 
and  the  amplitude  of  oscillation.  The  results  given  by  Garrick 
in  Reference  3  then  represent  the  performance  in  the  limit  of 
vanishingly  small  oscillations.  However,  the  effect  of  the  slip¬ 
stream  is  not  negligible  since  it  is  of  the  same  order  as  the 
thrust.  For  this  case  of  pure  heave  the  mean  thrust,  power,  and 
efficiency  are  given  by  the  following  equations. 
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[CL  ip  (00  c/2 ) 2  c 
a 

(ho*)U 

2  ( F2  +  G2) 

[ 10-AII ] 

[CL  ipO  c/2)3c 
a 

(ho*)U 

2  F/k 

[ 11 -All ] 

TUo  F2  +  G2 

f l-2k2 

(h  *)2 
v  0 

F-(F2  +  G2 )] 

J  [  12 -All ] 

?  '  p 

The  speed  at  which  these  occur  is,  of  course,  just  Uq  which  is 
given  by  Equation  [ 9— All ] .  These  quantities  are  plotted  in  Fig¬ 
ures  70,  71,  and  72  where  the  effect  of  the  slipstream  can  be 
seen  for  the  cases  of  hQ*  =  0,  0.1,  1.0,  and  10.0.  For  comparison 
the  results  obtained  from  the  large  amplitude  analysis  are  also 
shown . 
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